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We study, both analytically and numerically, the dynamics of elastic boundaries such as crack 
, fronts in fracture and surfaces of contact in solid on solid friction. The elastic waves in the solid give 

On ' r i se to kinks that move with a characteristic velocity along the boundary. As stopping kinks pass 

0^ ' through they cause moving parts of the boundary to stop. Starting kinks cause stationary parts of 

the boundary to move. We study the interaction of these kinks with disorder that arises from the 
spatial variations of the friction constant or fracture energy. In the absence of elastic waves, elastic 
boundaries with disorder operate at a critical point leading to a power-law distribution of slip events 
and self-affine boundaries. Elastic waves result in relevant perturbations at this fixed point. Slip 
\Q ' events beyond a critical size run away and the velocity of the boundary jumps to a nonzero value 

when the external load is increased above a threshold. We analyze in detail a class of simple models 
that capture the essential features of bulk vectorial elasticity and discuss the implications of our 
' results for friction and fracture dynamics. 

Moving elastic boundaries occur in a wide range of problems. Examples are the interface of contact in solid on 
i-j— | ' solid friction, the crack front in fracture, or the line of contact in the case of peeling of an elastic film off a substrate. 
, There remain several unsolved puzzles. For instance, why the frequency of the earthquakes averaged over all faults 
is proportional to a power of the seismic moment [Qj, while some faults show an excess of large events. Power law 
, t \ statistics of slip events may imply that the boundary is rough on all length scales just before it begins to move. 
d . Experiments on cracks restricted to move in a plane between two pressure welded Plexiglas sheets find that the crack 
front is self affine near onset of fracture [Q . The dynamics of an earthquake fault during a single slip event has received 
much attention since the inversion of seismic data by Heaton || discovered slip bands of finite spatial extent. Nielsen 
et al Q found corroborating numerical evidence in simulations of faults with rate and state dependent friction. Both 
studies use acoustic models, neglecting the vectorial elastodynamics of the bulk. Details of how these slip bands arise 
q ■ and how earthquakes start and stop are yet poorly understood. 

Analytic progress in problems with moving elastic boundaries is difficult and numerical simulations are intensive. 
Fisher and collaborators JB| considered a class of models to study rupture along heterogeneous earthquake faults. 
These were found to naturally operate at a critical point leading to power-law scaling of earthquake statistics. They 
also considered dynamical effects arising due to friction laws and elastic waves in an acoustic model. A variety of other 
' simplified models to study the dynamics of boundaries omit some important features such as long range elasticity || , 
bulk elastic waves and the vectorial nature of the problem |||| ■ To reproduce the complexity in the dynamics of 
~H ■ friction for instance, the transition between stick slip and steady sliding, Carlson and Batista introduced complicated 
\ friction laws ||. 

In this letter, we show that vectorial elastodynamics of the bulk results in complex friction dynamics even when 
the boundary obeys Coulomb friction. We find kink solutions that move along the boundary with a specific velocity 
and identify generic features of the equations of motion for elastic boundaries that would support such solutions. To 
understand the interactions of kinks with heterogeneities, we study a class of models, in the spirit of Re f. which 
capture the effects of elastodynamics. In our models, these kinks qualitatively alter the dynamics of the boundary 
""^J ■ near the onset of motion. 

To motivate our study, let us consider the problem of an elastic solid sliding on and always in contact with a 
flat infinitely rigid substrate. Let y be normal to the substrate and z be the direction of sliding. We study a two 
dimensional problem by requiring translational invariance along the x direction. Let f(z,t) be the displacement in 
the z direction of the solid's boundary relative to the substrate. An external shear stress <x^ = T and a normal load 
cr^ = P are applied far away. We obtain an effective equation of motion for f(z,t) by requiring that on the boundary, 



J c yz = /icr yy , when the solid is sliding, / > 0, 



(1) 

<7y Z < liOyy, when it is stuck, / = 0. 

(Note that the solid can slide only in the positive z direction, hence / > 0). We assume that the friction coefficient fi 
is independent of slip or sliding velocity Jul. Given f(z,t) and that elastic waves that radiate into the bulk are not 
reflected back, one can solve the equations of elastodynamics to obtain the resulting stresses at the boundary. Then, 
Eq. (|l|) becomes 

i^ZM =£(*,«) 0(£(*,t)), (2a) 
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where, 



E(s,t) 



dz'dt' J(z -z',t- t') df ^f ) + T - 7, 



' z',t'<t 

with 7 = ^P, and J = J lr i C tion = Jyz ~ M-/yy> where 



dt' 
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, for 1 < \z\/t < y/n, 



and 



Jyy " 2^ «9z 



Z /ty/z 2 /t 2 -ly/l-Z 2 /Kt 2 

l + (z 2 /t 2 -l)(l-z 2 /nt 2 ) 



for 1 < |z|/i < \//t. 



(4) 



Here k is the square of the ratio of the longitudinal to the transverse bulk wave speeds. A plot of the friction kernel 
■/friction as a function of time t at a fixed z is shown in Figs. [I] (c) and (d). 

When T > 7, there exists a unique steady Jll| state with the whole boundary moving at a constant velocity. We 
now seek, for T < 7, steady state solutions of Eq. (|J), f(z,t) = F(z — st) — F(£). Substitution into Eq. (|J) with 
J = ./friction yields 



A(s)F'(Z) + B(s)P 



(£-e) 5 



d£,' + T- 1 \Q[-sF\0\ =0, 



(5) 



where V denotes the Cauchy principal value and prime denotes derivative with respect to £. The coefficients in Eq. 
are given by 



A(s) 
B(s) 



s 2 sign(s)(l — s 2 1 K)\fs 2 



1 



2[l + ( S 2 -l)(l-s 2 / K )] 



s 2 /k \p sign(s) ys 2 



1 



27r[l + (s 2 -l)(l-s 2 /«)] 



(6) 
(7) 



We have shown p2| that a solution to Eq. (||) exists if and only if there exists sq such that -B(so) = 0. Steady state 
solutions can therefore propagate only when /i > fi c — \j\fk — 1 and at a velocity so = — Wl + l//i 2 - Furthermore, 
since so is negative, the expression in the curly brackets in Eq. (||) must vanish when F' > 0. For loads below 
threshold, i.e., T < 7, we find a steady state solution of the form 



This is a starting kink moving in the negative z direction. A stopping kink F'(!;) 

no 



- 7-T 



7-T 



A(s ) 
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[1 — O(^)], and pulses, 



A(so) — ®(£ — a )]' an y w idth a are also solutions. These predictions are confirmed by our numerics. 

These kink solutions are stable. When the speed of the backward moving kink \s\ > \sq\, B > and the integral 
term in Eq. (||) diverges logarithmically at the kink tip, £ = 0, so as to slow it down. If, on the other hand, \s\ < \sq\ 
the same term diverges with the opposite sign speeding up the kink. 

The motion of a planar crack through a three dimensional elastic solid is also described by Eq. (||) with J = 
Jf r acturc(z, t) sketched in Fig. |l| (d). Here, f(z,t) is the shape of the crack front. The equation of motion is the 
statement that the local elastic energy flux to the crack tip is spent to create new surfaces there. The driving force T 
is related to the external loads and 7 is the surface energy . Here again we find starting and stopping kinks that 
move with a characteristic velocity. Unlike in the case of friction, these kinks move in either direction along the crack 
front. 

Further progress can be made by realizing that the motion of elastic boundaries is generically described by Eq. (|^) . 
It possesses several remarkable features that arise due to the elastodynamics of the bulk and lead to the existence 
of steady state kinks. First, the equation for the motion of the boundary is first order in time. Second, the kernels 
J(z,t) have the following properties: (i) J(z, t < \z\/c) — 0, where c is the longitudinal wave velocity, (ii) the kernel 
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is homogeneous degree —2, so that J(z, t) = ^ Z J^ with j(0) = 1, and (iii) J is non-monotonic as a function of f at all 
z. In particular, at a fixed z, J is negative for some period after the arrival of the longitudinal sound wave and then 
changes sign . This change of sign implies there exists a speed sq < c such that for all z 

z|/s 

J(z,t')dt' = (9) 

\z\/c 

when J(z,t) = J(—z,t). The equation (^|) is a sufficient condition for steady kink solutions that move with speed so 
in either direction to exist.jj] To see that, we first note that condition (ii) implies that the steady state equation (||) 
is generally valid. Furthermore, we can rigorously show that condition (m implies that B(so) = 0. Thus steady state 
kink solutions exist for all such kernels. 

The detailed form of kernels that arise in problems of moving elastic boundaries depend on various parameters such 
as the microscopic time and length scales, and also on reflection of waves from free surfaces, etc. We would like to 
consider a simple model that captures effects of vectorial elastodynamics to understand the interaction of the kinks 
with disorder near the onset of motion. We therefore study a model kernel, shown in Fig. |l| (a), 

Jmodd = i[-6(t - \z\) + (2 + a)Q{at - \z\) ~ ae{c 2 t - \z\)}, (10) 
z z 

with 1 > c\ > c 2 . At a fixed position z this model kernel is zero until time t = \z\. In the time interval t € (|z|, |z|/ci) 
it is equal to — 1/z 2 . It then jumps to (a + l)/z 2 until time t = \z\/cn at which it settles down to its static value of 
1/z 2 . 

For T < 7, the equation of motion (^|), with the model kernel, admits starting and stopping kinks that move in 
both directions. Their velocity can be calculated using equation (||). It is 

f Cl(Q + 1) , for ( Q + l)( Cl -c 2 )>c 2 (l- Cl ) 

so = \ a + 2 - % (11) 

; — , otherwise. 

(2 + a — Cx — aci/c 2 

This is indeed confirmed by our numerics. 

Let us now turn to the effects of disorder which arises due to spatially varying material properties, such as the 
coefficient of friction or the surface energy. We introduce disorder by making 7 in Eq. (|J) position dependent. It is 
important to note that in general, the function 7 depends both on the coordinate z along the boundary as well as the 
position of the boundary /(z, t). In our simulations, the variable 7(2, /) is uniformly distributed in [0, T]. All lengths 
are measured in units of the correlation length of 7. 

We study Eq. (||) numerically by a finite difference method. Let us first consider the case the model kernel with 
a = 0. We increase the external driving force in small increments from zero. These increments are chosen in such a way 
as to dislodge the most weekly pinned point of the boundary. Avalanches begin with two starting kinks that nucleate 
there and move outward until they encounter a tough patch. Stopping kinks are then nucleated and move inward 
until the boundary is arrested stopping the avalanche. The dynamic exponent z is defined by the duration r ~ i z 
of an avalanche of size £. The avalanches show a power law distribution up to a characteristic length £_ with longer 
avalanches being much rarer. Below the length scale £_ the boundary is self affine with (\f(z, t) — f(z', t)\ 2 ) ~ \z — z'| 2< % 
which defines the roughness exponent, As the load approaches the critical load, T t qs , diverges as (T t qs — T)~" , 
defining the correlation exponent, v. Above this load, the boundary as a whole begins to move with an average velocity, 
v ~ (T — T t q ^)^, f3 being the velocity exponent. 

The Eq. (0) has been studied extensively for a quasistatic model where sound waves in the bulk are neglected, so 
that J = J qs (z,t) = 1/z 2 |]l5| . Note that J m odci(a = 0) < J qs for all z and t, and for t > I/C2, the two are equal. The 
no-passing rule p6| implies that the static exponents C an d v are the same in both cases. For the quasistatic model, 
the physics is determined by a critical fixed point and a 2 — e renormalization group calculation yields, with e = 1, 
C, = 1/3, v = 3/2. The dynamic exponents for our model arc different from the quasistatic model and a renormalization 
group calculation yields the dynamic exponent, z — 1 and (3 = 1. These results agree with our numerics, and the 
results are shown in Figs. [2] and ||. 

We find that for a > the structure of avalanches is the same as before for small enough loads. The typical 
avalanche size increases with the load. When it reaches a critical value £ c determined by the parameters a, c\ and C2, 



1 When the kernel is asymmetric, kinks may only move in one direction, as in the case of friction. 
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the starting kinks zip through the whole system causing it to move. We make this inference from the fact that the 
threshold load T t (a, ci, C2), is lower than that for the quasistatic kernel T t qs and that 

T t qs -T t (a, Cl ,c 2 ) ~a 3 / 2 ( Cl - c 2 )g( Cl ), (12) 

as shown in Fig. (|4|) . This scaling form of the suppression of the threshold load is in agreement with our analytical 
results [[l2|. Thus, as in the acoustic models in Refs. [p]]l7|, we expect a power law distribution of small events with 
the same power law as the quasistatic model combined with characteristic system size events. Also the velocity jumps 
to a finite value at T — T^(oe, c±, C2) as shown in Fig. |^. From our study we expect qualitatively similar behavior near 
threshold in the case of friction and fracture. 

In conclusion, we have demonstrated that vectorial elasticity of the bulk must be taken into account in order to 
properly describe dynamics of elastic boundaries. Specifically, steadily propagating kinks arise naturally in elasto- 
dynamic models. Pulse solutions propagating with a specific velocity, where parts of the boundary are moving can 
be constructed from starting and stopping kinks. Spatial width of the pulse solutions is arbitrary and must be set 
by some nonlinearity in the system or the reflection of waves from free boundaries. These pulse solutions should be 
observable in experiments on fracture similar to those of Schmittbuhl ||. These kinks can presumably be triggered 
at an edge of the sample. Simulations of solid on solid friction including the bulk elastodynamics jl8| suggest that 
the key features of the kernels obtained from a continuum treatment persist and hence kinks solutions must still exist 
though the velocity of motion of these kinks will be determined by the detailed shape of the kernel. Preliminary 
results indicate that the effects of sound damping in the bulk is to set a spatial extent for the pulses, but the details 
remain to be explored. More complicated friction laws, such as different static and dynamic friction coefficients, lead 
to complex behavior left for future studies. However, we expect kink solutions and the qualitative threshold behavior 
to persist. 

In systems with disorder, the dynamics near the onset of motion is controlled by nucleation of starting kinks by 
weak patches and of stopping kinks by tougher regions. Once an avalanche reaches a critical size which is set by the 
details of the kernel, it runs away spanning the entire system. Therefore, there should be an excess of system size 
avalanches beyond the crossover scale. Even when one includes the finite response time of a point on the boundary to 
changes in the stress, the same physics persists. Since avalanches beyond a critical size run away, rare weak patches 
should play an important role that has yet to be explored. 

In this letter we restricted our attention to one dimensional boundaries. When the spatial dimensionality of 
the boundary is two, " kink-on- kink" solutions arise. Thus, a possible mechanism for stopping a spreading circular 
earthquake, is a stopping kink-on-kink that propagates around the earthquake front. Questions about the dynamics 
of these kinks as well as their interaction with disorder remain open. However, models in which kinks are featured as 
the fundamental excitations offer significant advantages over full scale three dimensional numerical simulations and 
should be studied. 
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FIG. 2. Boundary statistics < f(k)f(—k) > just before it begins to move for the model kernel with no bump (a = 0). 
The solid line is the fc~ 5//3 scaling expected from the renormalization group arguments. Length are measured in terms of the 
correlation length of the randomness. 
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FIG. 4. Threshold force suppression for the model kernel with respect to the quasistatic value versus the 3/2 power of the 
bump height a. The forces T are in natural units. Diamonds correspond to C2 = 0.417, squares to C2 = 0.333, and circles to 
C2 = 0.2. The inset shows the threshold force as function of the ci — C2. In both graphs ci = 0.5 is fixed. 
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